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Abstract 



In récent years, the use of variational analysis techniques in lattice QCD has been demonstrated to be remarkably successful in the 
investigation of the rest-mass spectrum of many hadrons. However, due to parity-mixing, more care must be taken for investigations 
of boosted states to ensure that the projected corrélation functions provided by the variational analysis correspond to the same states 
^""^ at zéro momentum. In this Letter we demonstrate a novel method for ensuring the successful and consistent isolation of boosted 



states by expanding the operator basis used to construct the corrélation matrix to include operators of the opposite parity. 
baryon, excited state, variational analysis, non-zero momentum, lattice QCD 



'Keywords 



(yQ Introduction 

One of the greatest successes of lattice QCD has been its ap- 
' I / plication to hadron spectroscopy. The rest masses of not only 
the ground states, but also many excited states, can be extracted 
I iwith relative ease through a combination of effective mass tech- 
^-iniques and variational analysis. For example, there has been 
fH significant récent success by the CSSM Lattice Collaboration 
I ^C-Hl, the BGR Collaboration fS], and the HSC Collaboration 

' Once an understanding of the spectra is obtained, the logical 
çvq progression is to investigate the structure of thèse hadrons, and 
again lattice QCD has provided the tools needed for précise de- 
termination of quantities such as the Sachs form factors, charge 
radii, and magnetic moments. To date, however, such analyses 
have primarily focussed on the ground states and the current- 
induced parity transitions. 

. Key to lattice QCD's ability to investigate hadronic structure 
• • Is the computationof two- and three-point corrélation functions 
^ for each hadronic state of interest at both zéro and non-zero 
final state momenta. While the zéro momentum two-point case 
corresponds the rest mass analysis and is well understood, at 
. ^."non-zero momentum more care must be taken with baryons to 
ensure the correct corrélation functions are extracted, especially 
when investigating excited states. 

In this Letter, we investigate the use of variational analysis 
techniques to extract corrélation functions for excited states of 
baryons at non-zero momentum. In Sections |2] and |3] we de- 
scribe the traditional approaches and highlight how states of the 
opposite parity can intrude into the analysis. Section|4]demon- 
strates a novel method for overcoming this shortfall. This 
method will be central to future baryon form-factor calculations 
involving excited states, for example, electromagnetic structure 
and transition analyses. 



The results presented in Section |4] are calculated on the 
PACS-CS (2 + l)-flavourfull-QCD ensembles 10], made avail- 
able through the ILDG ISJ. They are 32^ x 64 lattices with 
P - L90, and employ an Iwasaki gauge action with non- 
perturbatively (9(fl)-improved Wilson quarks. In particular, we 
demonstrate proof of principle on the ensemble with the heav- 
iest quark mass, with k^ ^ = 0.137, corresponding to a pion 
mass of 702 MeV. 

2. Ground State Analysis 

To extract the mass for the ground state of a baryon, consider 
an interpolating operatoiQ^' that couples to this baryon and con- 
struct the corresponding two-point corrélation function. 



^(p;f) = 2e-'P''(Q|;fWJ(0)|Q). 



(1) 



By inserting a complète set of states I = 2 |B; p; s) {B; p\ s\ 
between operators and noting our use of Euclidean time, we 
can rewrite this as 

^^P;^) = ^^e-^«<p"(Q|;r(0)|B;/^;^)<B;p;^IHO)|îî), 

B s 

(2) 

where Eb{ç) - (ni^ + IpP)'^^ is the energy of state \B) at mo- 
mentum p, and niB is its mass. 

Let us now define coupling constants Ab and Àb such that 



(Q|^|B ;/?; s) = Àb* — UB*{p, s) and 



niB* 



<B+; p- s\x\0.) = Ab* J'^ubAp, s) 



(3) 



'Corresponding author. 



In gênerai thèse operators are not definite in parity, so that while they may 
transfomi as positive- or negative-parity operators, they couple to both positive- 
and negative-parity states. 
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for positive-parity states \B^}, and 



Mb- 



{n\x\B ; p; s) = Àg- — 75 ug-ip, s) and 



(B - p; sm = -Ab- y Ûb-(p, s) 75 (4) 

for negative-parity states \B^}, where mb± and mb± are positive- 
parity Dirac spinors. This substitution allows us to write the 
spin sum as 



^Mb±(/7, s) ub4p, s) 



y ■ p + niB* 
2mn± 



(5) 



and we can split Equation ^ into positive- and negative-parity 
components: 

At this point, we construct a new corrélation function by tak- 
ing the spinor trace of Q with some Dirac matrix F; 

G(r;p;f) = tr(r^(p;f)). (7) 
In particular, by taking F = ;= (7° + I)/2 we get 



B* 

+ ^ ÀB- ÀB- e 

B- 



_e^_,Eb- -niB- 



Er- 



, (8) 



and by taking F = F" := (7" - I)/2, 

G(r-,p;t)^yAB.AB.e-^-'^^^^ 
rr-* Er* 



+ ^ ÀB- ÀB- e 



-T .-e,-,Eb- +mB- 



(9) 



At zéro momentum, niB - Eb, and so F- act as parity pro- 
jectors, allowing us to separate out and individually examine 
the positive- and negative-parity contributions. Furthermore, in 
the limit as f —> 00, only the lowest lying state will remain, so 
that 

G(F^0;r)'"S^'2^B.^B.e-"'«î', (10) 

where \Bq) are the lowest-lying, positive- and negative-parity 
states. This allows us to easily extract the masses of the 
positive- and negative-parity ground states by defining the 
positive- and negative-parity effective masses of a corrélation 
function G as 

At non-zero momentum, mb + Eb, and so do not act 
as perfect parity projectors; contaminations from states of the 



opposite parity are introduced because the boost to non-zero 
momentum begins to mix the parities. This situation was in- 
vestigated in ||9[], where a projector of the form 



r"(p) := T 



niB, 



£b:(P) 



, TT 



(12) 



was introduced to remove a single contaminating state, the low- 
est State of the opposite parity. If this is the only nearby con- 
taminating state, this projector allows for the définition of an 
effective energy in the same fashion as the effective mass in 
Equation (fTTT i: 



G(F^(p);p;0 large, 

G(F±(p);p;f+l) 



Eb-M- (13) 



However, if there is more than one nearby state contaminating 
the corrélation function, this large-f approximation will again 
fail. 



3. Excited State Analysis 

The situation is potentially more challenging when using 
variational techniques. Here, we consider not just an indi- 
vidual interpolating operator;i' in isolation, as in Equation ([TJ, 
but the n X « corrélation matrix ^(p; f) formed from the cross- 
correlation functions between n différent operators i = 
CHEIII- Thatis, 



gijip; t) := Yj e"'" <"l^,W;?/0)|Q> . (14) 

X 

As before, we can parameterise the corrélation matrix in terms 
of coupling constants À" and À" , this time of each operator (/) 
to each state (a), to obtain 



,-£„t , y • P + ^B* 



^,(p;0 = Ç^r^f e- 



2Er 



(15) 



and we can again take the spinor trace with an arbitrary Dirac 
matrix F to get 

G,v(F;p;f) := tr(F^,/P;0) = ^ If e-^"' C(F; B; p), (16) 

B 

where 



C(F; B^;p) tr F 



.y ■ p + mB± 



(17) 

2£b*(p) / 

encapsulâtes the kinematics of the state 

We define a set of "perfect" operators { (ff ] that perfectly 
isolate individual energy eigenstates, so that 

(Q|0"|B^; p; s) = ô"\P^ z" u''(p, s) and 
V Eb" 

{BP; p; s\m = ^1^?' ^(P- s). (18) 
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Using the linearity of the operator space, thèse can be written 
as linear combinations of the original operatorsQ, with 



r-Y.u'^Xi and cf>"^Yj<X'- 



(19) 



1=1 



1=1 



Then, acting on Equation ( fTSI l from the right with the coefficient 
vector u" := {u", . . ., mJJ') and similarly from the left with : = 
(Vp . . . , v"), we obtain 



G(r; p; f) u'' = À" r e-'^" ' C{T; a\ p) and 



v" ' G{T; p; t) = z" /l" e"^" ' C{T; a; p) 



(20) 



where /l" :- {À" , . . . , À"), and similariy for A" . Moreover, act- 
ing on G(r; p; f) with both and v" at the same time pro- 
duce a diagonal, fully eigenstate-projected corrélation matrix, 
v""" G(r; p; f) = ô"l^ f ? s'^' ' C{T\ a; p), and hence we can 
construct corrélation functions that contain single energy eigen- 
states by 

^(np;?) := v"^G(r;p;f)u'' = z"7 s'^'' C{T;a;p). (21) 

With the t dependence in Equation (1201 constrained to the 
exponents, it is natural to write récurrence relations in / such 
that 



G(r; p; t + At) u" = e"^» ^' G(r; p; t) u" and 



v"^ G(r; p; f + Af) = e"^» ^' v"^ G(r; p; f) 



(22) 



That is, the coefficient vectors u" and v'' are the right and left 
generalised eigenvectors of G(r; p; t + At) and G(r; p; t), with 
generalised eigenvalue e"^" 

Importantly, for sufficiently large «, the lowest-lying states 
will preferentially be isolated, with ail higher excited state con- 
taminations constrained to the highest few states extracted. In 
extending an analysis from considering just the ground state, as 
in Section |2] the naïve approach is to continue to use the same 
Dirac matrix F in the parity projection, under the expectation 
that ail states isolated would be of a single parity. Indeed, at 
zéro momentum this is the case, since the parity projection is 
perfect. 

At non-zero momentum, however, the states of the opposite 
parity again contribute to the parity-projected corrélation mat- 
rix, and it is therefore possible the corrélation matrix analysis 
will isolate thèse contaminating states instead. The modified 
projectors of Equation (fT2t can remove only a single state of 
the wrong parity, leaving any other low-lying states of that par- 
ity alone to continue contaminating the analysis. Consequently, 
anoffier approach is needed to ensure the reliability of vari- 
ational analysis techniques in extracting individual energy ei- 
genstates. 



-This assumes that the set [Xi] of original operators spans the operator 
space. In reality, we are unlikely to have enough operators to completely span 
the space. In this case, the equality becomes an approximation, albeit improv- 
ing as the number n of operators used increases. The approximation can also be 
improved by perfonning the analysis at later times t and ï + Aïin Equation j22t 
since the Euclidean time évolution suppresses excited state contributions. 



4. Using an Expanded Operator Basis 

Instead of attempting to remove the contaminating states of 
ffie opposite parity, such as through the modified projectors, and 
analysing the positive- or negative-parity channels individually, 
we analyse both channels simultaneously using a single, ex- 
panded corrélation matrix. This expanded matrix is formed by 
recognising that iîxi is an interpolating operator that couples to 
the sector of interest, then so is -y s Xi (except it transforms with 
the opposite parity). Moreover, ffie sets {^i } and {ysXi I are 
pairwise linearly independent (indeed, at zéro momentum the 
subspaces spanned by each set are disjoint). Correspondingly, 
the adjoint sets are {xi ] and { -p, 75 ). 

If the Qijiç; t) of équation (fT4l l corresponds to the corrélation 
function formed from operators Xi and^j, then the correspond- 
ing corrélation function for y^Xi and^^ is 

G,j(9;t) ^Yj^-'""' {Q.\y,xii^)xm^) 
^^e-iP^rs {Q.\Xi{x)Xjm^) 
= 75 GuiV\ t), (23) 

where we have introduced a prime on the index to dénote the ex- 
tra 75 in the interpolator Similarly, the corrélation functions for 
Xi and -Xj 75, and y^Xi and -Xj 75, are g',/ (p; t) = -^//p; t) 75 
and ^,vj'(p; f) = -75 ^,j(p; f) 75, respectively. Then, using the 
invariance of the trace opération under cyclic permutations, for 
a Dirac matrix F the associated "projected" corrélation func- 
tions are, in the same order, 

G,v/F; p; t) = tr(F 75 Gijiç- 1)) = Gij{T 75 ; p; f), (24) 
Gif (F; p; = - tr(F ^,/p; t)y, ) = -G,v(75 T; p; f), (25) 



and 



G,v(F;p;r) = -tr(F75 ^//p; Ors) = -G,;(75r75;p;f). (26) 
Hence, we can write the expanded matrix in block form as 



G'(F;p;f) 



G(F;p;0 
G(-75r;p;f) 



G(F75;p;0 
G(-75F75;p;f) 



(27) 



This matrix will not be singular if the original « x « matrix 
is also non-singular, and hence we can use it to isolate twice 
as many states to compensate for the présence of the additional 
states of opposite parity at non-zero momentum. Since the op- 
erators should couple more strongly to states of the same parity 
as that which they transform as (in the case of F^, from Equa- 
tions (O and (|9]l the effective coupling goes as £ -H m ~ 2m for 
states of the same parity and E - m ~ for ffie opposite), this 
expanded corrélation matrix is nearly block diagonal. Hence, 
the eigenvectors corresponding to a state of one parity should 
be nearly orthogonal to those of states of the opposite parity, as 
expected. 

We can gain insight into the amount of "leakage" between 
parities by considering ffie éléments of the eigenvectors corres- 
ponding to operators of the opposite parity. In Figure [T] we plot 
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# 


D (1 U ) 


loi (GeV) 


1 


(0,0,0) 


0.000 


2 


(1,0,0) 


0.427 


Q 

J 


^^ 1 n\ 
U , i , (J^ 




4 


(1,1,1) 


0.740 


5 


(2,0,0) 


0.854 


6 


(2,1,0) 


0.955 


7 


(2,1,1) 


1.046 


8 


(2,2,0) 


1.208 


9 


(3,0,0) 


1.281 



Table 1: Momenta used in this analysis. 



the eigenvector associated with the positive-parity, ground state 
A as a function of momentum-squared, and in Figure |2] the 
eigenvector associated with the negative-parity A(1405). Here, 
we use the 12x12 expanded corrélation matrix formed from the 
operators ™d x'^^ (Equations (2.5a), (2.5b), and (2.7) 

of 11211 ') with 16 and 100 sweeps of gauge-invariant Gaussian 
smearing Hl, where T = r+ (I + /)/2. 

We can clearly recognise the effect of Lorentz contraction in 
increasing the overlap of less smeared operators at large mo- 
menta, as demonstrated in lll4ll . In both Figures [T] and |2] the 
dominant two same-parity components correspond to the same 
flavour structure — for the A and;^f^-' for the A(1405) — ^but 
with différent amounts of smearing. At low momentum, the 
dominant smearing is the 100 sweeps, while as the momentum 
increases the contribution from the smaller 16 sweep source 
increases. We suspect the two discontinuities in this trend, 
between the fourth and fifth, and eighth and ninth momenta, are 
due to the lack of SO(3) rotational invariance on the lattice; they 
correspond to the transitions from momenta diagonal across the 
lattice to momenta purely in the ;ic-direction. The momenta used 
are tabulated in Table [T] 



5. Conclusion 

While the relative sizes of the eigenvector components asso- 
ciated with opposite-parity operators are small (< 5%) com- 
pared to the dominant same-parity component, they are non- 
etheless non-zero and, if neglected by only using the traditional 
non-extended corrélation matrix, prohibit the clean isolation of 
energy eigenstates. As we have demonstrated in this Letter, by 
extending the corrélation matrix to include operators of the op- 
posite parity, we can investigate and cleanly isolate both parity 
channels at the same time. Moveover, since the parity trans- 
formation amounts to simply pre- or post-multiplying (or both) 
the projection matrix F by yg, it is a computationally inexpens- 
ive procédure. 
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Figure 1 : (colour online) Momentum-squared dependence of the eigenvector associated with the positive-parity, ground state A. 
Components associated with operators of the same parity are indicated in the left plot while components associated with operators 
of the opposite parity are illustrated in the right plot. At zéro momentum, ail opposite-parity components are consistent with zéro 
within errors. Eigenvectors are normalised so that they have unit 2-norm and positive maximal-absolute-value component. 
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Figure 2: (colour online) Momentum-squared dependence of the eigenvector associated with the negative-parity A(1405). Com- 
ponents associated with operators of the same parity are indicated in the left plot while components associated with operators of the 
opposite parity are illustrated in the right plot. At zéro momentum, ail opposite-parity components are consistent with zéro within 
errors. Eigenvectors are normalised so that they have unit 2-norm and positive maximal-absolute-value component. 
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